Reentrant Phase Diagram of Two-Dimensional Electron Gas 
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We study the phase diagram of an electron on a square lat- 
tice in a transverse magnetic field. Using a renormahzation 
scheme, we show that the inequality of the two next-nearest 
neighbor couplings destroys the fat critical regime above the 
bicritical line and replaces it with another reentrant extended 
phase. Furthermore, the universal strange attractor describ- 
ing self-similar fluctuations of the localized phase is replaced 
by a new fixed point. The reentrant extended phase and the 
new fixed point of the localized phase belong to the univer- 
sality class of an anisotropic triangular lattice. 
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As the properties of the Harper equation ||], which 
describes an electron on a square lattice in transverse 
magnetic field, are fairly well understood much less 
is known about more general quasiperiodic tight binding 
models (TBM). In this paper, we study the generalization 
which results from taking into account both the nearest- 
neighbor (NN) and next-nearest-neighbor (NNN) inter- 
action in the electron problem. The associated TBM has 
the form ||] 

ta(V'/c+l + V'fc-l) + 2tf,COs[27r(fccr + 

+ exp[z2^(A:CT + (f))]{tabe'^^ ^k+i + V,e-^""Vfe~i} 



+ exp[-i27r(fcCT + (t>)]{tabe''"'ipk-i + t^ie' 



= E^k (1) 



Here ta and tb are the NN couplings while tat and 
are the diagonal NNN couplings. The quasiperiodicity 
results from the irrational parameter a describing the 
magnetic flux. The Harper equation corresponds to van- 
ishing of the NNN couplings. Furthermore, in the limit 
t^b equal to zero, the model 

Recently, the above model was studied in the isotropic 
limit where tab = tab- Although the wave function 

tpk in general is a complex function of the lattice index 
k, Tpk can be taken to be real in the case. These studies 
revealed a very interesting phase diagram : For 2tab < ta, 
the model belonged to the universality class of the Harper 
model with both extended (E) {tb < ta) and localized (L) 
{tb > ta) phases and a critical (C) point at ta — tb where 



the system had the full square symmetry. On the other 
hand, for 2tab > ta the model was found to belong to a 
new universality class where there was no E phase but 
instead the C phase existed in a finite parameter interval 
tb < 2iafc. For tb > 2tab the states were exponentially 
localized. The E and the C phases were separated by a 
bicritical line 2tab — ta- 

Detailed decimation studies |^ showed that the wave 
functions within the fat C phase above the bicritical line 
were self-similar (at the band edges) only at certain spe- 
cial values of the parameters. These special points corre- 
sponded to universal limit cycles of the renormahzation. 
However, for generic parameter values, the fractal char- 
acteristics of the critical wave functions did not exhibit 
self-similarity and were conjectured to be described by a 
strange attractor of the renormahzation flow. 

In another interesting study of this model, it was shown 
that the fluctuations of the wave functions in the L 
phase mimiced the behavior in the C phase The L 
phase of the Harper universality class was described by a 
renormahzation fixed point of the strong coupling limit 
tb/ta 00 while the L phase above the bicritical line 
was described by a strange attractor of the associated 
renormahzation 1^. 

The general case of the model where tab and tab are 
not equal and the resulting TBM is complex has not 
been fully investigated. In the three parameter space, 
A = tb/ta, a = 2tab/ta, and (3 = 2tab/ta, using the 
duality property of the model, Han et al. Q calculated 
the Lyapunov exponent of the model analytically and 
concluded that the system is localized for A > 1 if 
{a + (3)/2 < 1 and for A > (a + f3)/2 otherwise. Apart 
from the existence of metal-insulator transition, noth- 
ing is known about the scaling properties of the complex 
model. 

In this paper, we study the complex model using our 
recently developed decimation scheme. We will confine 
ourselves to the case where a — (V5 — l)/2. It is shown 
that the phase diagram changes discontinuosly as a — /3 
becomes different from zero. The C phase of the real 
TBM above the bicritical line is replaced by another E 
phase. This complex reentrant E phase is described by 
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a strong NNN coupling of the triangular lattice and is 
shown to be related to the weak coupling limit of the 
Harper equation. Furthermore, the invariant strange set 
of the renormalization describing fractal characteristics 
of the fluctuations in the L phase degenerates to a fixed 
point of the renormalization. These fluctuations are de- 
fined by the equation p| 



-l\k\ 



(2) 



where 7 is the Lyapunov exponent which vanishes in the 
E and C phases and is positive in the L phase. In other 
words, rjk is equivalent to the original wave function ipk in 
the E and C phases whereas in the L phase rjk describes 
the fluctuations around the exponentially decaying wave 
function. Knowing the analytic formula for the Lyapunov 
exponent, it is easy to write a TBM for rjk, resembling 
Eq. (1), in the L phase 

We will use a decimation approach to describe the scal- 
ing properties of rjk, the wave function in the E and C 
phase and the fluctuations of the wave functions in the L 
phase, for a speciflc value of energy. In our studies below, 
we will focuss on the quantum state with minimum en- 
ergy Emin- In addition to flxing the quantum state, one 
has to also fix the phase factor </) to a critical value in 
Eq. (1) so that the wave function remains finite asymp- 
totically. The nondivergent wave functions are needed 
to determine the scaling properties as has been discussed 
previously. 

The key idea of the decimation scheme is to connect 
the wave function rj^ at an arbitrary site k with two 
neighboring Fibonacci sites k + -Fn-i-i and k + Fn where 
F„+i=Fn + Fn-i |,|: 

{k)v{k). (3) 

The way how the decimation functions /„ and e„ are 
placed in the decimation equation is somewhat arbitrary 
but here we choose the form which causes the asymptotic 
limits of the decimation functions e„ and /„ as n — > 00 to 
be bounded in all three phases. The additive property of 
the Fibonacci numbers provides exact recursion relations 
for the decimation functions e„ and /„ (|:^,^: 



e„+i(fc) = 
fn+i{k) = 



l + A/„(fc) 

fn-l{k + Fn)fn{k + F,,) 



l + AUik) 

A = e„_i(fc + i^„) + /„_i(fc + F„)e„(fc + F„) 



(4) 
(5) 



For fixed fc, the above coupled equations for the deci- 
mation functions define a RG flow which asymptotically 
{n 00) converges on an attractor. The C phase is dis- 
tinguished from E phase by the existence of nontrivial 
limiting behavior. With anisotropic NNN coupling, the 
attractor is a p-cycle in all three phases for E — Emin 



]lOt . The asymptotic p-cycle for e„(0) and /„(0) deter- 
mines the universal scaling ratios 



(j = lim r]{Fpn+j)/r]{0); j = 0, 1. 



(6) 



whose absolute values are equal to unity in the E phase 
and less than unity in the C phase [p^ . 

In Fig. 1 we show the phase diagram obtained by an- 
alyzing the asymptotic behavior of the decimation func- 
tions in different parts of the parameter space. As soon 
as a and /? differ, a — /3 is an irrelevant parameter and the 
phase diagram is determined solely by two parameters: 
A and {a + /3)/2. In the parameter range (a + (3)/2 < 1, 
the decimation functions become asymptotically real ap- 
proaching the same universal cycles as for the Harper 
equation a = /? = . For (a -I- /3)/2 > 1, a 7^ /3, the 
decimation functions stay complex also asymptotically. 
The cycle length p in this case is six in all three phases. 
However, considering the absolute value of the decima- 
tion functions (or the scaling ratio C), one observes a 
3-cycle on the line {a + /3)/2 = 1 and a flxed point for 
{a + f3)/2 > 1. The 3-cycle observed at the critical point 
C is different from the universal 3-cycle observed along 
the rest of the (bi)critical line AC. Moreover, these flxed 
points are all universal and do not depend on the actual 
parameter values. Table I summarizes various universal- 
ity classes of the model. It is interesting to note that 
the line CL divides the localized phase into two different 
universality classes. That is, the scaling properties of the 
self-similar fluctuations in the region BCL are different 
from those of the region LCE. Furthermore, the bound- 
ary line CL has its scaling characteristics different from 
the two regions that it separates. 

The reentrant E phase described by a complex 6-cycle 
with complex scaling ratios is not as trivial as the real 
fixed point of the E phase below the bicritical line de- 
scribing the weak coupling limit of the Harper model. 
The latter depends neither on the phase 4> nor on the 
lattice index k and can be easily solved from a fixed 
point equation: /„(fc) = a and e„(fc) = —a^ resulting 
in ^ = 1. In the complex E phase, complications arise 
from the fact that the decimation functions do depend 
both on (j) and k. But the absolute value of a decima- 
tion function has the same constant value as in the real 
case. Noting the fact that the decimation functions and 
the scaling ratios above the line ACL remain the same, 
the reentrant E phase can be understood in the limit of 
highly anisotropic triangular lattice a — > 00, P/a — s- 0, 
and A/a ^ 0. In this hmit, Eq. (1) reduces to the TBM 
describing the weak coupling limit of the Harper model. 



Cfe+i + Ck-i — - — Cfe, 

tab 



where Ck is related to ipk via 



Ck = exp(i27r0fc) exp(i7r(Tfc^)'i/'fe. 



(7) 



(8) 
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Hence, the wave function of the infinitely anisotropic tri- 
angular lattice is related by a phase factor to the ex- 
tended wave function of the similar NN square lattice. 
The above equation shows that the scaling factors Q 
exist only for special (i.e. rational or those related to 
the golden mean) values of the phase (f>. The relation 
Fn(T ~ Fn-i — (— cr)" impHes that 



O'F^ — Fn-lFn 



i-iy 



1 + 2(7 



(9) 



From this and Eq. (8) it follows that for (/) = 0, 1/2 
Cj = ±exp[±i7r/(l + 2a)] 



(10) 

The above scaling is relevant also in other parts of the 
complex E phase because the RG flow is attracted by 
the same 6-cycle in the whole region. Taking advantage 
of the above limiting solution, it is possible to derive an 
explicit expression for the 6-cycle. Substituting Eq. (8) 
into the decimation equation (3), we obtain 

e„(fc) = eJi(fc) exp(-i27r(/)F„) exp[-i7ra(i^2 ^ 2fcF„)] 
/„(fc) = /„^(fc)exp(i2^0F„_i) 

X exp[i^(7(F„\i - i^,? + 2fci^„_i)], (11) 

where ejj and are the decimation functions corre- 
sponding to the TBM (7) describing the weak coupling 
limit of the Harper equation. From these equations we 
see that e„ and /„ are functions of the fractional part of 
fc(T, denoted by < fccr >, only. Therefore, we can write 
the decimation functions in terms of the renormalized 
variable x = (— tr)"" < /ccr > ||^. For simplicity, let us 
assume that 4> ~ 0,1/2. Applying the relation (9) and 
the fact that 



(-ct)"F„ 



(-1)" - 0-2" 



(12) 



l + 2a 

we obtain six different limiting function pairs of the form 

e*(x) = ±0-2 exp[±i7r(2j; - 1)/(1 + 2cr)] 
f*{x) = ±crexp[±i7r2(cra; + 1)/(1 + 2a)] (13) 

as n tends to infinity. These pairs form a 6-cycle of the 
recursion (4-5), written in the continuos variable x 

f \ Aen{-ax) 
'--'^^^^- l + AU-ax) (14) 
r t N /n-i(g^a: + a)fr,{-ax - 1) 

^"+^(")^ l + AU-ax) (1^) 

A = en-i{a^x + a) + fn-i{a^x + a)en{-ax - 1). 

For the above 6-cycle, 1 -|- Af*{—ax) = a. 

Therefore, the characteristic feature of the reentrant E 
phase is the fact that the decimation functions depend 
explicitly on x. These functions are complex and con- 
sequently the universal scaling ratio has both real and 



imaginary parts but the absolute value of ( is unity. This 
is unlike the Harper E phase where the real universal 
functions are site independent and are given in terms of 
the powers of the golden mean. 

In summary, we have shown that with anisotropic NNN 
couplings, the renormalization behavior for the TBM 
describing an electron on a square lattice is a lot sim- 
pler than in the isotropic case. Firstly, the renormaliza- 
tion strange set corresponding to the fat C phase in the 
isotropic case is replaced by an attracting cycle associated 
with trivial scaling properties (i.e. E phase). Secondly, 
the bicritical lines of the isotropic case remain critical also 
when the NNN couplings are not equal but the renormal- 
ization attractor is again simpler (i.e. a cycle). Thirdly, 
the fluctuations of the exponentially localized wave func- 
tions are described by a universal fixed point and not by 
an infinite strange set as in the isotropic case. 

The novel feature of our model is the existence of two 
extended phases separated by a critical line signaling a 
transition from real scaling ratios to complex ones. The 
two extended phases respectively fall into the universal- 
ity classes of the weak coupling limit of the NN square 
lattice and the strong NNN coupling limit of the trian- 
gular lattice. It is interesting that also on the other side 
of the localization boarder (i.e. in the L phase), the scal- 
ing behavior is divided into three different universality 
classes. The Harper L phase corresponding to the strong 
NN coupling is separated from the strong NNN coupling 
phase by the CL line with its own scaling properties. 
Therefore, for a fixed value of NNN coupling, as the NN 
coupling A is varied, the behavior beyond the localization 
transition appears to be shadowed by the behavior before 
the localization transition: the existence of E phase be- 
fore localization results in a nontrivial fixed point in the 
L phase while the existence of a period p limit cycle de- 
scribing self-similar critical states before the localization 
onset leads to the appearence of period p limit cycle of the 
renormalization after localization describing self-similar 
fluctuations. 

The phase diagram discussed here describes the univer- 
sal properties of square and triangular lattices . There- 
fore, we believe that our results can be experimently re- 
alized on two-dimensional mescoscopic systems. |L1| 
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FIG. 1. The reentrant phase diagram of the TBM (1) 
for the case a ^ (5. The (bi)critical line AC separates the 
two extended phases while the critical line CE describes the 
onset to localization. The extension of the AC line (line CL) 
also separates the localized phase into two different regions 
with different scaling properties. For the corresponding phase 
diagram in the isotropic case a = (3, see Fig. 1 of ref. [5]. 

TABLE I. The universal scaling ratios for various univer- 
sality classes at the band edge. In order to avoid showing all 
p scaling ratios, some of which differ only by signs of the real 
or imaginary part, we show all different (|7?e(Cj)|, \Im{Q)\). 
(jic = 1/2 for all other parts of the parameter space except the 
line AC where (j)c varies as a function of A. In the reentrant 
E region ACE, the absolute value of C, is unity. Note that the 
last three lines show the scaling properties of the fluctuations 
in an exponentially localized wave function. 
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